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ABSTRACT 


The capacity of an air terminal for Short Takeoff and Landing air- 
craft is analyzed. The terminal is considered to be operating as part of 
an intra-urban air rapid transit system. The air traffic flow through 
the terminal is modeled by a computer simulation written in both the 
FORTRAN IV and GPSS languages. The model is used to solve the 
traffic capacity problem under two sets of traffic control rules. In the 
first case, existing FAA rules, which require 3 miles separation between 
oP ivals and 2 miles between an arrival and a departure, are used. Ina 
second case, the rules are 2 miles between arrivals and 1 mile between 
an arrival and adeparture. A detailed description of the modci is 


presented so that others might use the model. 
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Lew INTRODUCHION 


One of the many problems facing our cities today is the lack of effec- 
tive and economically feasible rapid transit systems. Solutions to this 
problem cannot be thought of in terms of today's hardware and technology 
because of the long lead time required for development of complex rapid 
transit systems. In fact, today one must consider systems for operation 
five or ten years in the future. 

One system that is currently being examined is based on the use of 
Short eo Lt and Landing (STOL) aircraft in this very short haul intra- 
urban market. This type of transportation system may be suitable for an 
urban area suchas the Sana meisco Day area, It womwld eomre Weems 
rather than replace existing systems, suchas San Francisco's BART. 
To the extent possible, the STOL system would use existing facilities. 
Where no air terminals now exist, suchas in downtown San Francisco, 
terminals would be built especially for the STOL system. 

ime System is envisioned to operate more like a bus line than an air 
line. There would be no cabin attendants and stops at terminals would 
be a matter of minutes instead of hours. Only one basic type of aircraft 
would be used and the system would be operated under a local monopoly, 
probably by a municipally owned company. 

The National Aeronautics and Space Administration has initiated 


both in-house and contract studies of the STOL rapid transit system. 








NASA examined the economics of the system [1], and the system was 
found to be economically feasible. Reassvured that the idea had merit, 
NASA let contracts for a detailed system design to the Boeing corpora= 
tion. The Boeing study proposed an aircraft to be used in the system 
and a terminal design and placement scheme. [2] The STOL port proposed 
had a single runway with multiple passenger gates. It was to be used 
only by the STOL system with no other traffic allowed. 

Boeing performed an economic analysis of their system proposal 
and found it to be economically feasible under the proper assumptions. 
The critical economic areas were the fixed costs of construction and 
operation, and the costs of unproductive ground and holding pattern time 
of the aircraft. [2] All of these costs are directly related since a greater 
amount of unproductive aircraft time requires more aircraft to satisfy a 
given demand than would be required if the unproductive time were short. 
The higher number of aircraft would require more terminal facilities to 
handle them. The cost of having to obtain and operate the increased 
number of terminals and aircraft could make the system economically 
unsound. Thus it is important to have both aircraft and terminals 
Operating efficiently to minimize the number of each required to operate 


the system. 








Il, STATEMENT OF THe Peer e mM 


The matter of facility capacity seems to be an important key to the 
potential success of the system. This paper presents the development 
of a model which can be used to determine STOL terminal traffic capacity 
in the context of the rapid transit system described above. 

Given a particular set of assumptions, a measure of capacity will 
be defined and capacity determined under two sets of traffic control 
rules. The two solutions will be made to demonstrate the effect on 
Eeeacity of ailowing reduced aircraft separation. 


mairst the problem will be solved using existing FAA traffic separa] 
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anarrival and a departure. The problem will then be solved using 
separation requirements of 2 miles between consecutive arrivals and 

1 mile between an arrival and a departure. This latter set of separation 
requirements reflect rule changes that may result from the use of more 


precise navigation and air traffic control equipment. 








ils, SOLUTION TECH Wews 


The economic success of the STOL rapid transit system depends 
Iteavily on the rapid and efficient cycling of airezaft throupm thertermimais. 
lime seonsequences of an inaccurate capacity figure therefore maken 
educated guess, or trial and error solution, unacceptable, since long 
aircraft delays caused by overloaded terminals or idle facilities caused 
by under utilized assets would prove disasterous to the system budget. 

Computer simulation was chosen as the solution method in light of 
the Mee escity of a solution and the obstacles posed by direct analytical 
P-caniques. The interdependencies between em aE events in the system 
Mould make direct solution an impossible task. If the dcpendcencics and 
distributions involved were simplified to the point that direct mathematical 
solution was possible, the model resulting would probably be only a gross 
mepresentation of the system. 

The simulation model described in the following sections was used 
in the problem solution. Five hours of simulated time were simulated to 
Peemmne the period covering evening rush hour, this being the time of 
heaviest traffic. The model was run for two hours of simulated time at 
ny 2 peak flow. The traffic flow rate was then increased over a period of 
90 minutes to peak flow and held there for one hour of simulated tine. 

The model then simulated another 30 minutes ata traffic flow of 1/3 


peak. 








Data was examined for the final 3 hours of the 5 simulated hours. 
The peak flow was adjusted until the traffic capacity was found. The 
criterion used to define capacity was the expected delay time for those 
aircraft which had to wait to start their approach. This was estimated 
using the average holding pattern time. When the expected delay time 
for those aircraft that had to wait reached 3 minutes, the airport was 
considered to be at capacity. This admittedly somewhat arbitrary 
criterion was chosen because it represents the author's opinion of the 
iaaexiinum amount of time a traffic controller could be expected to be 
able to delay an arrival by indirect routing - hence it is the point at 


which an arrival would have to enter a holding pattern. 
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[VV Aso UMP FrOnNs 


Pee LIST OF ASSUMPTIONS 

The assumptions that were made in the solution of the problem were: 

1. The aircraft involved is a multi-engine STOL aircraft with fast 
load/unload capability. 

2. The aerodynamic stall speed is 60 MPH in the landing 
configuration. 

er 6Line Speed the aircrait flies on the approach is normally distripe 
Wiced Path mean 82.3 MPH and standard deviation of 5.04 MPH. 

4, The landing speed is normally distributed with mean 77 MPH 
and standard deviation 4.¢5 MrT. 

bemeeelhe aircraft decelerates after landing at a constant rate ot 
10 Wises’ umtil a speed of 10 MPH is reached. At that time, theater. 
turns off the runway and is clear in 5 seconds. 

6. The distribution of time between arrivals is exponential and all 
arrivals are independent. 

7. Taxi time between the runway and the gates is 1.5 minutes for 
arrivals and departures. 

8. Gate time is exponentially distributed with a mean of 3 minutes. 

we «6cLhe time required to position thevaireratt in position tor take— 


off on the runway is .25 minutes. 


Jill 
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10. The take-off 1s modeled as a constant accelerations nema 
standing start at the rate of 10 Rees untila speed of (7 (ewe er 
Beached. 

li. The terminal has an approach course that is 9 miles long. 

12. 13 miles must be flown back to the start of the approach course 
if a go-around (or wave-off) is necessary. 

more Lhe aircraft is capable of executing a safe wave-oil tromaean 
point on the approach course greater than 0.1 miles from the landing 
ean. 

14, The terminal is equipped with 3 passenger gates. 


IS. The weather is good with no wind. 


Pee oLiFICATION Cr ASSUMPTICNS 

This section is included to explain the reasons for the assumptions 
made inthe problem solution. 

1. The aircraft was modeled after the proposal made by the Boeing 
Semporation. [2] 

2. The speed distributions for both the approach speed and landing 
speed were based ona study of speed distributions of present day pas- 
senger jets [3], and the premise that the STOL aircraft will be able to 
be flown as consistently near designed approach speed as present day 


Memiuarercial aircraft. There ig reason to believe they will in fact be 


easier to fly than today's jets. 
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3. The inter arrival time assumption is based on studies of actual 
arrivals at commercial airports [4,5,6]. One may doubt the validity of 
drawing an analogy between contemporary commercial airports, with 
many independent (in the statistical sense) airlines, and the one company 
pmol port. However a parallel situation does exist. Instead of ditierent 
airlines, we have aircraft on different routes inthe system. If one 
considers each route as a different line, the exact same traffic situation 
exists at the STOL port and a commercial air port. 

4. The gate time distribution assumption is based on the characteris- 


tics of similar types of short stop systems [4]. 
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VV. SENSITIVITY ANALYesis 


This section is included to indicate the sensitivity of the solution to 
the assumptions made. 

The speed distribution assumption is conservative. The previous 
section indicated that the normality assumption would probably result 
in speeds more variable than would actually exist. This is stated to be 
conservative because a smaller variation in speeds than that assumed 
would result ina larger peak flow at capacity. This happens because 
the PEeitic separation is attained at the start of the approach. If two 
@omeccutive arrivals have different approach speeds and the faster one 


be 1 


fomlows the slower one, the Second may be tco close pening the 


Lees L “& 
edt ee . 


Ch 
0) 


teeaue. if this happens, the first aircraft will not be clear of the run- 
way and the second will have to go around. 

The solution is entirely dependent on the mean approach speed 
assumed. Since the separation required is specified as distance between 
aircraft (rather than time between aircraft), aircraft with faster approach 
speeds will be able to attain this distance in less time than slower air- 
craft. The wind plays a role in this also, since the ground speed 
metally determines the time needed to get the required spacing. An 
indication of the effect of ground speed can be seen by noting that the 
time needed to get the approach separation is a linear function of the 


ground speed. At the assumed ground speed of 82.3 MPH, a capacity of 
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22 aircraft per hour can be handled using the 3 mile separation rule. 
Applying the ratio of peak flow to speed, one can predict a peak capacity 
of 24 aircraft per hour at a speed of 94 MPH, and 19 per houriet 2a 

The selection of the landing speed and deceleration assumption 
Pecome critical only when the time required for an airerait tovcleanuine 
runway aiter landing is so great thai there is no opportunity to inserta 
departure between two arrivals. Naturally the amount of time required 
for a departure is as important as the time needed fora landing. The 
time needed for a departure to taxi into position for takeoff is not critical 
since this can be done while the previous arrival is decelerating. What 
is required then is that the time between arrivals be such that a de- 
celeration from landing speed to taxi speed and an acceleration to take- 
off speed can occur. Looking at the mean speeds involved, we find that 
a landing followed by a takeoff takes about 25 seconds. The average time 
between arrivals (using 2 miles separation and 82 MPH) is 88 seconds. 
There is obviously a lot of slack time available, which means that the 
Bepoiwne and takeoff procedures are not critical to the flow rate. 

Taxi time has no effect on the waiting time and simply adds on to 
Piestotal cycle time. Thus as far as the capacity of the terminal is 
@emeerned, the taxi time is really immaterial. 

iae pate time could have an effect om the solution. With the 3 
minute mean time assumed, the gate queue was normally empty for all 
flow rates tested. Obviously, the amount of gate time available (number 


of gates multiplied by the time period being considered) can be such that 


1S 








the gates become a bottleneck, the solution then would have to consider 
waiting time at the gates as well as at the holding pattern. 

The weather assumption is not important. The model operates air-~ 
craft under Instrument Flight Rules. These procedures will operate the 
same regardless of the actual weather, unless the weather is so bad it 
forces closure of the terminal. It is true that bad weather results in 
higher landing speeds and lower deceleration rates, but since there is so 
much slack available in the landing and roll out phase, these longer run- 


way times will have no effect on the solution obtained. 
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Vi. ANALYSIS OF Daw 


As has been stated, the statistic used to indicate capacity was 
average conditional waiting time. If this time was 3 minutes or greater 
the airport was considered over capacity. The model operates on the 
basis of abstract time units, with the interpretation in this solution 
being that 1 time unit equals 5 seconds. Thus the question was really 
to find the maximum flow rate that resulted in an expected conditional 
waiting time of 36 time senaslte or less. As apreliminary step, based on 
ine eeotral Limit Theorum, the assumption was made that the average 
conditional waiting time ve distributed normally with mean equal to 


s 9 
Va Tiance UNnKNnoOWwrTi. 


Ci, 


the true mean of the waiting time an 

With the assumptions stated above, the following procedure was 
used to determine the maximum peak flow rate that resulted in a mean 
conditional waiting time of 36 time units or less. 

ie OA value of peak flow was picked and the model run for 10 
replications. 

2. a. If the average value of the mean conditional waiting time 
was greater than 36, a test of hypothesis using the 95 percentile point 
ofthe '"t'' distribution was made. Ho was that the true mean was less 
than 36, Ha that the true mean was greater than or equal to 36. 

b. If the average of the mean conditional waiting time was 


less than 36, the test of hypothesis was made with Ho being that the 


Rey, 








true mean was greater than 36, and Ha that the true mean was less than 
or equal to 36. 

3. a. If Ho was not rejected, more data was taken andapooled 
with the previous data. The test of hypothesis procedure described above 
was repeated. 

b. If Ho was rejected, the peak flow rate was adjusted. If the 
average of the sample tested was greater than 36, the peak was decreased. 
If the average of the sample just tested was less than 36, the peak flow 
rate was increased. The model was run for 10 replications at the new 
peak flow rate and the above procedure repeated starting at step 2. 

The flow chart of Figure 1 illustrates the procedure just described. 
This search and test procedure was used to try to minimize the computer 
time needed to solve the problem. Since the actual flow rate capaci 
was not known within 10 increments, testing all reasonable possibilities, 
using a sample size large enough to insure a desired confidence level, 
would have been an expensive procedure. As it was, approximately 45 
minutes of computer time was required to get the data that was used. 

A disadvantage of this procedure is that, due to the sequential nature 
of the data collection and hypothesis testing, it cannot be said that the 
hypotheses were tested at the 95% level, even though the 95 percentile 
point of the ''t'’ distribution was used to define the critical regions. 
Lindgren [7] presents a discussion of the problems associated with 
sequential sampling. The true value of the type I error in this case 


(the probability of rejecting a true Ho) is greater than 5%. Unfortunately, 
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if is probably beyond calculation. The difficulty of caleulaticongiceec ass 
pounded by the fact that the sequential tests were not of the same sample 
size. Often a different number of data points resulted from different 
batches of runs. This was a phenomenon peculiar to the monitor system 
of the computer installation where the model was run. If the run time 


exceeded a certain real time limit, the run was terminated at that point. 
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Vil. RESULTS 


The model was run under the conditions as stated. For each run of 
5 hours of simulated time, the following was recorded: 

i Lhe proportion of aircraft that did not have to waite staan 
Pmeir aporoach. 

2. The average waiting time for those aircraft that had to wait to 
start their approach. 

bee che average time for an aircraft to cycle through the termina 
Peasured from arrival at the start of the appreoaciro completiense: 
farce Olt, 


Tht =< 7 ee Cf Aathe- RieaG 
Tie three scts of data Wer 


a) 
x) 
al) 
a 
(iP 
jos 
i$ 
> 
") 
2 
2] 
) 
+, 
> 
DD 
en 


Semrelation between them. The criterion decided upon to measure 
eapacimy, that is, the expected conditional waiting time, proved highly 
variable. Unfortunately, no precise correlation could be found which 
would allow the use of another measure to indicate that the desired 
delay limit of 3 minutes average delay had been reached. 


Tables 1 and 2 contain a summary of the results obtained. 


fad 
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Hg 


H1 


PS St 
SUMMARY OF RESULTS 


3 MILES BETWEEN ARRIVALS, 
2 MILES BETWEEN ARRIVALS AND DEPARTURES 


N PNW M S C He H1 
10 .28 54.10 20.21 202.76 R 
1] 33 45.20 19. 93 197. 14 R 
99 . 38 39.31 21.75 187.54 R 

105 43 Bee? | Sage 179.74 R 
48 46 31.88 ere 178.24 R 


imeak arrival rate imeaitcratt per nour. 

Number of replications on which figures are bascd. 
Average proportion of aircraft that did not have to wait. 
Average conditional waiting time, in model time units. 
The sample standard deviation of the conditional waiting time, 
Mmastmode! tile units, 

fincwaverave Cyclic Mine, 10 Ie@dc meme mimalic. 

An "R'' in this column indicates that the hypothesis ''the true 
mean is greater than 36'"' was rejected. 

An '"R" in this column indicates that the hypothesis "the true 


mean is less than 36'' was rejected. 


ZZ 





TABLE 2 
SUMMARY OF RESULTS 


@¢ MILES BETWEEN ARRIVALS, 
1 MILE BETWEEN ARRIVALS AND DEPARTURES 


R N PNW M S C Hg H1 
6 Tr . 30 318 023 elec 192. 48 R 
34 167 35 34.96 45 AD 187.28 

oem 123 36 32.54 19.50 185.30 R 

32 «50 . 42 26.14 14,99 178.87 R 


Column headings in this table are the same as those used in Table I. 


al 


* The flow rate of 34 per hour could not be resolved. That is, the 
hypothesis that the true mean was greater than 36 could not be rejected, 
even with the large number of data points obtained. In view of the large 
sample variance at this level, and the closeness of the sample mean to 
36, it was decided to abandon this flow rate and be satisfied with knowing 
that the flow rate just below it was under capacity, while the one above 
it was over capacity. 

The flow rate of 35 per hour is omitted because it could not be 


represented by an integer number of time units of 5 seconds width. 
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VIII. CONCLUSIONS 


It was concluded that under existing FAA traffic separation rules, 
the capacity of the terminal was a peak flow rate of 22 aircraft per hour. 
Under the revised separation requirements of 2 miles between arrivals 
and 1 mile between an arrival and a departure, a peak flow of 33 aircraft 
merenour can be handled. 

There is a significant advantage to be gained from allowing the 
closer spacing. It would certainly be worth while to carefully examine 
the separation rules to be applied. If is toovexpensive in fermicvoislocs 
Sapemty to require unnecessary distance between aircraft. 

Bie capaci igWres resulting from thit solution are considerably 
more pessimistic than those used in the Boeing study. Their assumptions 
resulted in a much more disciplined arrival rule than the Poisson process 
assumed here. In fact, they assumed that the aircraft arrived at the 
Initial Approach Fix already properly spaced for the approach [2]. As 
a result, they were able to sustain a flow rate of 28 aircraft per hour 


under existing rules, and a flow rate of 41 aircraft per hour under the 


2 mile - 1 mile rules. 
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IX. MODEL DESCEI Eien 


The model was designed to represent the specific scenario ofa 
STOL terminal operating in the environment of an intra-urban rapid 
transit system. Asa result, the following assumptions are built into 
the model and cannot be changed without re-writing the program. 

ime Lhe model considered traffic to be hoimogenious elm taiscueocee 
homogenious applies both to the aircraft type involved and to the method 
Sasoperation. All aircraft are considered under positive radar control 
and eeroting under imstrument flight rules (17 R)- 

2. The terminal modeled has a single runway with a single approach 
eourse from the hoiding paticern to the runway. 

fe Lhe holding pattern is at the initial approach fix (JAF): 

a lleaweues are Of the tirst-1n=t)cusoutlyoce 

ee Arrivals are given priority Over departures. 

6. Aircraft that have had to go around on an approach are given 
priority for another approach over other aircraft in the holding pattern. 

mihe traffic flow in the model is as indicated in Figure cd. Figure > 


shows the physical system that is modeled. 
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XA. THE GPSS PR@Gee 


The simulation was written using General Purpose Simulation 
System (GPSS). This system generates entities and moves them through 
a sequence of queues and facilities according to user defined rules. It 


uses abstract time units to measure flow rates and delay times. 


fee GPSS ANALOGIES 

The analogies that are used in the GPSS model are as follows: 

1. Queue 1 represents the holding pattern, queue 3 the gate queue, 
and queue 4 the takeoff queue. 

mee Hacility | represents the approach course, facility 2 the 7ume 
way, and facility 21 the runway threshold. 

3. wotorage 1 represents the passenger gates, with the storage 
capacity determining the number of gates to be simulated. 


mee each GPSS entity represents an aircraft in the system, 


Peer oo FUNCTIONS 

lmeorder to use the model, one must have a basic knowledse of 
So functions since they are used to describe transit and delay time 
probability distributions. To define a function, the user must input a 
set of ordered pairs to the GPSS program. The program interprets the 
first number of each ordered pair as the independent variable, and the 
second as the dependent variable. The entire function is approximated 


by linear interpolation between the poin’s described. 
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The GPSS program generates a uniform random variate on the 
interval 0 to 1 (or 0 to 1000) to 3 significant figures. By using this 
number as the independent variable, a function value is obtained from 
the user defined function. This value is tnen multiplied by the mean of 
am@eecdistribution in question to pet a realization trom the distribuivom 
lawouner words, the function defined ieenteente a mapping from a uniform 
distribution to the desired cumulative probability distribution, nor- 
malized so that its mean is l. 


The exact mechanics and form for inputting the ordered pairs will 


mewcescribed in section All. 


fees OUTPUT 

See Goo Sutpu inized inko sections as follows? (he weades 
icetind it helpful to refer to the GPSS output from the sania in 
meppendix B). 

1. Source listing of the GPSS deck. In this listing each operational 
card is assigned a block number and each card in the program is assigned 
Meeara number by the compiler. It will be noted that these program 
assigned numbers bear no relation to the card numbers punched in 
columns 78-80 of the source deck. 

2. Apartially compiled source listing (this listing is omitted in 
mppendix B). 

3. The flow surnmary. This is the first of the statistical output 


data from the run. This summary is a count of the total number of 
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entities that have moved through each block, as well as the number 

that werepresent in the block at the end of the run. Note that if statistics 
are kept from other than the start of the run, the total shown for each 
block is summed only over the time for which statistics were kept. 

4, The runway utilization statistics. Three values are shown 
here. The average utilization is computed as the proportion of time 
the runway was occupied. The number of entries is a count of the 
number of operations, takeoffs and landings, that occurred on the run- 
ieee ihe average time per transaction is the average time each run- 
way operation took, in GPSS time units. 

5. The runway threshold utilization statistics. The data in this 
section is computed in the same manner as the runway utilization 
statistics. 

6. The gate statistics. The capacity simply indicates the number 
of gates that were simulated. The average contents shows the average 
number of aircraft that were at the gates. The average utilization is 
the proportion of available gate time that was used. The average time 
per transaction is computed as the average time an aircraft remained 
at the gate, in GPSS time units. The current contents and maximum 
contents are self explanatory. 

7. The queue statistics. These are self explanatory except for 
the zero entity concept. A zero entity is one which did not have to 
wait in the queue, but simply passed through. The $AVERAGE TIME/ 
TRANS column shows the average time each entity had to wait in the 


queue, not counting zero entities. 
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8. Table of transit times. This 1s 2 tabulation cistmes. oO. epies 
each aircraft took to transit the model, in GPSS time units. The table 
structure, that is the range and increment size, is user defined. The 
overflow portion of the table includes all values greater than the table 
upper bound. All intervals in the table have as their lower bound the 
upper bound of the previous interval. The first interval has zero as 


its lower bound. 
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XI. THE FORTRAN PROGRAM 


The FORTRAN section of the model is used solely as a pre- 
processor to convert input data toa form acceptable to the GPSS 
program. It accepts inputs in a convenient form and units, suchas 
MPH, miles and Mince. It outputs printed output for the information 


Mimtee user, and punched output to be included in the GP5Ss deck. 


PeeenpORTRAN INPUTS 

All inputs to the FORTRAN program are on punched cards. Appendix 
B shows a sample input deck. Inputs fall into two types, required and 
optional. Naturally, all required inputs must be assigned values for 
each run. The normal FORTRAN convention relating to variable types 
has been followed, all variables are real (decimal) numbers except 
those whose names start with the letters I, J, K, L, M, N, or O. 
Input variables which are real numbers must have a decimal point 
supplied, even if the value assigned is a whole number. Integer var- 
iables must not have decimal points included. 

mie first character of each card of the input deck must be punched 
Starting in column Z, while the last character on each card must be 
metore column 72. The first characters on card 1! of the input deck 
must be an ampersand followed by the word INPUT, followed by at 


least one blank. 
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A variable is assigned a value by punching the variable name, an 
equal sign, its value and a comma. An example is: ARRT=5.3, . After 
the last comma on the last card, an ampersand followed by the word 
END must be punched. 

As many cards as are necessary may be used. Spacing between 
variables is unimportant. The last character on all cards except the 
last must be a comma. That is, a variable definition must not be split 
between two cards. 

fia vector variable is involved, the entire vector may be input uy 
naming it and then listing as many values as there are vector elements. 
An example of a vector variable is NFLAG which contains 6 elements 
and is an integer variable. It may be inout by punching NFLAG=3, 3, 3. 
4,1,0, which sets NFLAG(1) to 3, NFLAG(Z) to 3, etc. The same 
result may be obtained by punching NFLAG(1)=3, NFLAG(2Z)=3, etc., 
erm BAG=3*3,4,1,0; where "'3*3" is ecuivalentito 3, 3,3: 

i Required inputs 

The following named variables are required inputs to the 
FORTRAN program. 
AoEaeD Distance in miles from the initial approach fix to 
the landing touchdown point. 
Je lee) aa), The mean approach speed of the aircraft, in MPH. 
ARRT The initial mean time interval, in minutes, 


between arrivals at the holding pattern. 
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AVGATE 


DS A 


Ue Pr ink 


LOWER 


LINC 


NGA TiS 


RUN 


Syeda Oh 


SED 


TIMEU 


WDIST 


The mean time required to load and unload 
passengers, measured sn anubes- 

The number of miles of separation required 
between consecutive aircraft on the approach course. 
The upper limit for the table of transit times, units 
ar Gatlinie limits. 

The lower limit for the table of transit times in 
time units. 

The number of time units in each tabular interval 
in the table of transit times. 

The number of passenger gates available in the 

ereae) eobiroretll 

The amount of time, in time units, to be simulated. 
Statistics are kept for this time period unless a 
value for STDY is specified (see Optional Input 
Ssectron tor ST DY): 

The number of miles of separation that is required 
between an arrival and a departure. 

The mean aircraft landing speed, in MPH. 

The number of seconds to be represented by each 
oe bese, 

The number of miles that must be traveled ona go- 
around to re-enter the approach course at the 


initial approach fix. 
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WOFFD The minimum distance, in miles from the touch- 
down point that a wave off can be safely initiated. 
2. Optional Inputs 
The following named variables are optional inputs to the FORTRAN 
program. The variabies for which a default value is listed will be set to 
that default value if the user does not assigna value to them. 
PoC Ch, The acceleration rate on takeoff, measured in 
ft/sec“. The default value is 10 tui coe 
NE LAG A six element vector used to flag eptiomsmn specs 
ifying probability distributions. Three values are 
meaningiul; !, 2, and 22 tne detauit valilerisme 
This variable is used to set up options ior the mol 
lowing variables: ARRT, APPSPD, TDSPD, AVGATE, 
TAXIIN, and TAXIOU. The preceding list is ordered, 
i.e. NFLAG (1) refers to ARRT, NFLAG (2) to 
APPSPD, etc. Each time the simulation needs a 
value to use for one of these parameters (the actual 
inter arrival time between two aircraft, for example) 
it realizes a value froma probability distribution 
according to the following rules: if NFLAG as 1 
(or default) a single value is realized with prob- 
ability one. This value is the mean value, as input. 
If NFLAG is 2, a uniform distribution is sampled. 
If NF LAG is 3, a user specified distribution is 


sampled. 
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AMODIF 


PACCEL 


CHANGE 


A six element vector used with NFLAG. If the 


corresponding NFLAG value is 1, AMODIF is 
ignored. If the corresponding NFLAG value is 2, 
the parent variable value (ARRT, APPSPD, TDSPD, 
etc.) is considered the mean, and the AMODIF value 
is half the range of the uniform distribution. That 
is, the distribution range is the mean + AMODIF. 

If the corresponding NFLAG value is 3, the AMODIF 
value is the number of points the user wishes to 
specify to define the function the simulation will use. 
The user then must supply the required number of 
function points to the GPSS program. To realizea 
value, the number assigned to the parent variable 
will be mmltiphed by the functicen valucucbsa 7a 

The program normally computes the time to accel- 
erate to takeoff speed on the basis of a constant 
acceleration at rate ACCEL from a standing start 

to landing speed. If this model is not satisfactory 
to the user, the number of time units required to 


accelerate can be input here. 


The option exists to change the mean interarrival 
time periodically during the simulation. If this 
option is to be used, the number of time units 
between changes in mean interarrival time is input 
here. Note that RUN divided by CHANGE should 


result ina whole number. 
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DECEL 


aig 


NPUNCH 


NWRITE 


IPP LAG 


ei) Y 


The deceleration rate after landing, measured in 
ft/sec”. The default value is 10 Cee 

This is a flag used to indicate that inclement 
weather is to be simulated. A value of 1 indicates 
this optiomis to be used. The result is the inereac= 
ing of landing 2-01] ovt bya factor of 1.1, 99 zene 
indicates noninclement weather is to be simulatec. 
The default value is zero. 

The logical value which is assigned to the unit to 
process punched output. The default value is 7, the 
normal IBM 360/67 punch code. 

The logical value which is assigned to the unit to 
process printed output. The default value is 6, the 
normal IBM 360/67 printer code. 

This isa flagvenairerait transiamtimes eA | Specities 
that each aircraft transit time is to be printed and 
tabulated. A zero specifies tabulation only. The 
default value is zero. This option should be used 
with caution. Each aircraft transit time printed out 
will résuit in three@ames of output. Tliis "ean easiiy 
lead to a large volume of paper. 

The option exists to keep statistics from other than 
the start of the run. Assigning a number to this 


variable specifies the number of time units from 





the beginning of the run at which statistics will be 
started. 

TAXIIN ihre mean number of miinwes required to taximtcem: 
the runway to the passenger gate. The default is 
2.0 minutes. 

TAXIOU Ihe mean numer of minutes required fo tasieiten 
the passenger gates to the end of the runway. The 
default value is 2.0 minutes. 

DRG oF The number of time units necessary to decelerate 
LOuta <a Speecadtter landino. | Uneadetaultcenaautes 
the time using a model of constant deceleration from 
landing speed to TURNOF (see below) ata rate of 
DECEL. 

UI RIN EF The speed in MPH at which the aircraft can tuemeotr 
the runway after landing and deceleration. The 
default value is 10 MPH. 

WINDD The wind direction relative to the runway heading, 
iWweasumed in decrees froma zero to) .o0, | ine detauls 
Values Zero, 


WINDS ihe wind “speed int viele [he defamievaiie 1s 2ere. 


me PORTRAN VARIABLES 
The following list contains the principle FORTRAN variables and the 


mares and quantities used to compute them. 
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Variable formats are required for some outputs. The vectors whicn 


Pre used to contain these formats are: FMT, BLANK, SFMT, DIGIT, 


erelo, STG, TAB, TABR, FPNCTS and 4! through Z15- 


HWIND 


CoP ED 


COMIT 


ATIME 


IHX (1) 


IXH(2) 


IXH(3) 


IXH(4) 


The computed headwind component. 

ihe mean ground speed. 1t is computedsby supe 
tracting the headwind component from the mean 
approach speed. 

The time before an arrival lands after which no 
more departures Gam pe cleared aheadvor hima... 
is computed by dividing the minimum allowable 
Spacing between arrivals and departures by the 
mean ground speed. 

The mean fime required for an approach. It is 
obtained by dividing the approach distance by the 
mean ground speed. 

The initial mean interarrival time converted to 
time units. 

Thevmean time needes@s#to obtain the required imvem= 
arrival spacing on the approach course. 

The mean time required to travel from the holding 
Paveemmerorrie coterie point. The cOmamit point is 
determined from the COMIT computation above. 
Ene maean time required to travel irem the commu 


point to the runway threshold. 
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IXH(5) 


IXH(6) 


ET | 


IXH(8) 


IXH(9) 


IXH(10) 


134s GA) 


IXH(12) 


1X (13) 


IX H(14) 


Perit 15) 


IXH(16) 


CHANGE 


‘The mean time needed to travel from the threshold 


to the touchdown point. 

The time required tovclear the runway after landing: 
It is computed by adding 5 seconds to the time 
required to decelerate airecr landing. 

THe meanstiime crequigce to taxi from thesnumnya, ale 
the gate. 

The mean gate time. 

WTiMie watean time meeded: totax) fostacmr nwa nei 
Pee oamcs 

ive stitiite sac ece cen tomtas] si homdeinc 8G @ istic tase 
way to position on the runway for takeoff. 

ihe timee required) toracecverare amemtaiceor 

The time required to travel back to the start of the 
approach if a go-around is required. 

The time from which statistics are to be kept ale 

GCL eit hae OMautimMicese or, 

A flag to indicate if variable mean interarrival time 
1s tonpe wsied, 

UnewmMeaneinec rarrivalesimerchange litera) ites 
Setee ater ne le INiGuiEe 

The number of gates to be simulated. 

The time intervals at which changes in interarrival 


time are to be made. 


4] 


IC The number of times the inecan inferanrive atime 
to be changed, if this option is used. It is computed 
by dividing RUN by CHANGE and subtracting 1 from 
the product. 

ithe program uses two subroutines called “FITIT” and "RFUNCiHS | 


These do format fitting and no other computations. 


fee ORTRAN OUTPUT 

Appendix B contains a sample of the FORTRAN output. The printed 
output consists of two types, informational and directive. 

The informational output gives the user a record of what the program 
has done. The first output produced is a copy of the values of all input 
Variabies. Some optional inputs which were not assigned values by the 
user may show negative values. This is normal and should not concern 
fre user. The program will print a copy of all punched cards produced. 

The program will check all 6 NFLAG values and their corresponding 
PBMODIF values. if an NFLAG is specified 2 or 3, and no AMODIF value 
is assigned, the NFLAG will be setto 1. All NFLAG values after the 
check is made will be printed. 

The directive output consists of instructions to supply data to the 
moo program. The data required will be a specified number of function 
points. All information is to be on punched cards. The program will 
specify the placement of the punched cards as "following card ‘ 


where the blank will contain a card number. The card numbers referred 


weeare punched in columns 77-80 of the GPSS deck. 
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XII. PROCEDURES FOR GPSS RUN 


This section contains the procedures the user should follow between 
the FORTRAN run and the GPSS run, as wellas some general remarks 
about the programs. 

A set of function points will be required for each user supplied 
mrobability function to be used (each NFLAG set to 3), and for mean inter- 
arrival time changes. For probability functions, the first number of the 
Sadered pair specifying the function point is the independent variable. 
This aan be realized from a uniform distribution. The range of the 
uniform distribution will be 0 to 1 for arrival rate, landing speed, gate 
Sefenetaxi in timc, and taxi oul time. It will be @ tc 1000 for approach 
meeea. Lhe second number of the ordered pair is the value of the depend 
Smeevatiable at the function point which is being described. 

mon inean interarrival time changes, the first number of each pair 
is an integer from 1 to IC (RUN/CHANGE -1). The second number is the 
mean interarrival time for that time period. 

Appendix B contains some examples of function definition cards 
miceluded inthe GPSS deck. All of these cards are punched starting in 
column 1 and not past column 72, with no imbedded blanks. As many 
cards as are necessary may be used. Ordered pairs are separated by 
the character /, and the numbers in each pair are separated by a comma. 


The last pair on each card is not followed by the character /. Decimal 








points may be omitted for whole numbers. The pairs must be ordered 
among themselves so that the independent variable is monotone increasing. 
Smmexample would be; 1,2.31/1.1,5/1. 11, 1054/2, 0.03. 

ie FORTRAN propram will) produce wumehediontouucy ic aeeronee 
included in the basic GPSS source deck. The cards will have numbers 
punched in columns 77-80 which indicate their relative position in the 
GPSS deck. Ifa like numbered card already exists in the basic GPSS 
deck, it is to be replaced by the new card. Otherwise, the new card is 
to be placed in the deck in numerical order. Gaps can be expected in 
the card numbering system. The GPSS deck should be restored to its 
@ricinal basic configuration as listed in Appendix A, before another run 
Hm Set up. 

The GPSS program will produce the same sequence of random 
muiImbers for each run. To obtain a different sequence of random num- 
bers, a card punched starting in column 8 with the word RMULT, and 
punched starting in column 19 with 2 commas and a number, may be 
included in the GPSS deck. The number must be a seven digit or smaller, 
odd number, and will be the seed value for the GPSS random number 
Benmeraror, This card should be placed immediately following the 
vid LATE card at the head of the GPSS deck. A different sequence of 
random numbers will be produced for each different RMULT number 


used. 
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The following pages show the input and output for a sample run of 
eae model. 

The first page shows the input cards that were used. Next is the 
FORTRAN printed output followed by the punched output. 

The GPSS source output shows the function definition cards that 
were included as per the instructions in the FORTRAN printed output. 
mae GPSS output has been edited to exclude the partially cornpiled source 
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